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ABSTRACT 

Transiting planet discoveries have yielded a plethora of information regarding the internal structure 
and atmospheres of extra-solar planets. These discoveries have been restricted to the low-periastron 
distance regime due to the bias inherent in the geometric transit probability. Monitoring known radial 
velocity planets at predicted transit times is a proven method of detecting transits, and presents 
an avenue through which to explore the mass-radius relationship of exoplanets in new regions of 
period/periastron space. Here we describe transit window calculations for known radial velocity 
planets, techniques for refining their transit ephemerides, target selection criteria, and observational 
methods for obtaining maximum coverage of transit windows. These methods are currently being 
implemented by the Transit Ephemeris Refinement and Monitoring Survey (TERMS). 

Subject headings: planetary systems - techniques: photometric - techniques: radial velocities 



1. INTRODUCTION 

Planet formation theories thus far extract much of 
their information from the known transiting exoplan- 
ets which are largely in the short-period regime. This 
is because transit surveys which have provided the 
bulk of the transiting planet discoveries, su ch as Super- 
WAS P (jPoUacco et, alJl2006fl and HATNet ijBakos et all 
2002), are biased towards this region. The two plan- 
ets which contribute to the sample of intermediate 
to long-period trans iting exoplanets are: HD 17156b 



(iBarbieri et all |2007[) and HD 80606b (jLaughlin et all 
I2009t iMoutou et all (2009h . the latter of which exhibits 
both a primary transit and secondary eclipse. In both of 
these cases, the detection was largely due to the inflated 
transit/eclipse probability caused by their extreme or- 
bital eccentricities. Both of these were observed to tran- 
sit through predictions based upon their radial velocity 
data. 

Planetary orbits may be considered in three basic cat- 
egories: short-period (< 10 days) planets, intermediate 
to long-period high-eccentricity (> 0.1) planets, and in- 
termediate to long-period low-eccentricity planets. The 
first type of planet are the focus of current studies, and 
we are beginning to gain insight into the second type 
(HD 80606b for example). Exploration into the struc- 
ture of those planets occupying the second and third orbit 
types will require probing parameter-space beyond that 
currently encompassed by the known transiting exoplan- 
ets. The science objectives of such an exploration include 
understanding how planetary properties, such as aver- 
age planet density, vary with periastron distance as well 
as providing the first observational data for exoplanct 
models with low incident s tellar r a diatio n. R ecent obser- 
vation s of HD 80606b by iGillonl ()2009f ) and iPont et all 
(2009) suggest a spin-orbit misalignment caused by a 



Electronic address: skane@ipac.caltech.edu 

1 NASA Exoplanet Science Institute, Caltech, MS 100-22, 770 
South Wilson Avenue, Pasadena, CA 91125 

2 Department of Astronomy and Astrophysics, Pennsylvania 
State University, 525 Davey Laboratory, University Park, PA 
16802, USA 

3 UCO/Lick Observatory, University of California, Santa Cruz, 
CA 95064 




n n-u-i-rl, r jT — n 



Periastron Distance [AU] 



Binsize = 0.1 AU 



Periastron Distance [AU] 

Fig. 1. — Histogram of periastron distances of known exoplanets. 
The shaded region corresponds to known transiting planets, the 
population of which is comprised of planets in short-period orbits 
and intermediate-period eccentric orbits. 



Kozai mechanism; a suggestion which could be investi- 
gated in terms of period and eccentricity dependencies if 
more long-period transiting planet were known. 

Figure [T] shows the distribution of periastron distances 
for the known exoplanets using data from Jean Schnei- 
der's Extra-solar Planets Encyclopaedia 4 . Also shown 
as a shaded histogram is the same distribution for the 
known transiting exoplanets. Clearly the periastron dis- 
tribution of the known transiting planets does not accu- 
rately represent the distribution of the entire sample of 
known exoplanets. Thus far, our picture of exoplanets 
is based on the planets with super-heated atmospheres 
through cither short-period orbits or intermediate-period 
eccentric orbits. Figure 1 shows how our view of plane- 
tary properties is highly biased towards planets of short 
periastron distance, thus demonstrating the need for de- 
tecting transits of exoplanets at larger distances in order 
to study atmospheres with greatly reduced flux from the 
parent star. 

Long-period transits will give us valuable insight into 
the structure of exoplanets that are more similar to 

4 http://exoplanet.eu/ 
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those in our own solar system. Additionally, host stars 
of long-period planets discovered using the radial ve- 
locity technique tend to be very bright, and poten- 
tially one of these could be the brightest star with 
a transiting Jupiter-mass planet. The brightness of 
such host stars can faciliate atmospheric studies via 
transmission spectra (iBurrows et alT 2006; Redficl d et al.1 
l2008t [Snellen et "all 120081) and thermal phase variations 
( Knutson et al.1 I2007D . The long period typically also 
translates to a longer transit duration, allowing signif- 
icantly more signal to be collected when trying to probe 
the atmospheric composition of these objects with high- 
resolution spectrographs. Since most bright late-F,G,K 
type stars have already been searched for short-period 
Jupiter-mass planets, the long period hosts are the obvi- 
ous sample to now search since their relative brightness 
and longer transit times are both very helpful in attempt- 
ing to probe the exoplanetary atmosphere. 

It is clear that the discovery of long-period transit- 
ing planets is very important. They arc, however, dif- 
ficult to detect in ground-based transit sur veys due to 
the t he presence of correlated (red) noise (|Pont et al.1 
2006) and effects of the o bservational window func- 
tion (|von Braun et al.l I2009D . Though combining data 
from various on-going transit surveys to search for 
tran siting long-period p lanets will improve detectabil- 
ity ([Fleming et al.ll2008l ). it is easier to discover tran- 
sits among the large number of long-period exoplan- 
ets already known from radial velocity surveys. Based 
upon Monte-Carlo simulations of the transit probabil- 
ities, it is expected that several of the known long- 
period planets should transit their parent stars due to 
the enhanced probabilities produced by eccentric orbits 
([Barnes! 120071 : iKane fc von Braunl l2008f l. There have 
been suggestions regarding the strategy for photomet- 
ric follow-up of these radial velocity pla nets at predicted 
times o f primary transit (iKanel l2007h and secondary 
eclipse ([Kane fc von Braunl I2009D . and t he instruments 
which could be used for such surveys (|L6pez-Moralesl 
2006). Here we describe a methodology through which to 
search for transits of known planets. In Section 2 we cal- 
culate transit windows for a large selection of the known 
exoplanets and show the impact of additional radial ve- 
locity measurements on refining the transit ephemeridcs. 
Section 3 describes the techniques through which op- 
timal target selection and observations of the transit 
window can be achieved. The methods described here 
have been successfully tested and implemented by the 
Transit Ephemeris Refinement and Monitoring Survey 
(TERMS). 

2. TRANSIT WINDOWS 

The transit window as described here is defined as a 
specific time period during which a complete transit (in- 
cluding ingress and egress) could occur for a specified 
planet. The limiting factor for successfully observing a 
known exoplanet host star during the predicted transit 
window is often the precision of the transit ephemeris. 
The quality of the transit ephemeris is primarily deter- 
mined by (a) the uncertainties associated with the fitted 
orbital parameters, and (b) the time elapsed since the 
most recent radial velocity data was acquired. Acquiring 
new high precision radial velocity data can easily miti- 
gate both effects. With a prompt photometric observing 



strategy after the orbital parameters have been revised, 
one can maximize the chances of being able to obtain 
complete coverage of the transit observing window and 
thus either confirm or rule out the transiting nature of 
the planet. 

2.1. Ephemeris Calculation 

The orbital parameters measured from fitting the ra- 
dial velocity data of a planet are sufficient for calculating 
a transit ephemeris. The predicted time of mid-transit 
can be calculated by using Kepler's equations. Firstly, 
the eccentric anomaly is calculated from the following 



E = 2 tan" 1 




(1) 



where e is the orbital e ccentr i city a nd / is the true 
anomaly. As described in lKand (|2007[ L the time of tran- 
sit mid-point will occur when lu + / = n/2, where u> is 
the argument of periastron. Substituting this for the true 
anomaly in Equation Q] thus yields the eccentric anomaly 
at the point of predicted transit. 

The mean anomaly, M, which defines the time since 
last periapsis in units of radians, is then computed by 

M = E - e sin E (2) 

which can be converted to regular time units using 

PM 

t M =— (3) 

where P is the orbital period. The predicted mid-point 
of primary transit can then be calculated using 

PM 



2vr 



nP 



(4) 



where t p is the time of periastron passage and the term of 
n x P incorporates the number of complete orbits which 
have transpired since t p . 

The uncertainties in the orbital parameters (assuming 
they are symmetrical) can be propagated through these 
equations to determine the uncertainty in the predicted 
transit mid-point, <5t m id, and the size of the transit win- 
dow, t wln . Note that these uncertainties will be equiva- 
lent to la if the orbital parameter uncertainties are also 
1<T. The size of a transit window is mostly dependent 
upon the uncertainty in the period and the time elapsed 
since last observations were acquired. Thus, the begin- 
ning and end of a transit window are calculated by sub- 
tracting and adding (respectively) the uncertainties in t p 
and P with respect to the transit mid-point, taking into 
account the number of orbits since periastron passage 
and the transit duration. The beginning of a particular 
transit window can be approximated by 

tbegin - (t P - St p ) + (P-8P)^-+ n(P -SP)-% (5) 

Z7T 2 

where 5t p and SP are the uncertainties in t p and P re- 
spectively, and td is the transit duration. Conversely, the 
end of the transit window is approximated by 

Wd = (tp + St p ) + (P + SP)^-+ n{P + SP) + ^. (6) 

Hence, the size of a given transit window is defined by 
subtracting Equation [5] from Equation [6l resulting in 



/wi„ = 2 I Stp + 5P^ + nSP j - 



(7) 
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Fig. 2. — Ephemeris calculations for the sample of 245 exoplanets for the first predicted transit after t p (see Section |2,2| |, These show the 
dependence of transit window (top-left), transit mid-point uncertainty (top-right), and predicted transit duration (botton-lcft) on period, 
as well as the relation between the transit window and the transit probability (bottom-right). 



which reduces to simply the transit duration as the un- 
certainties in P and t p approach zero, as expected. It 
is clear from Equation [7| that 5P has the potential to 
rapidly dominate the size of the transit window if the 
number of orbits since discovery becomes sufficiently 
large. Equation [7] may be re-expressed as follows 



SP = 



ir(t v 



td — 25t p 



M + 2nn 



(8) 



which can be used to determine the period uncertainty 
needed in order to achieve a certain transit window for a 
fixed 5t p and td (see Section [2~3l for examples). 

These equations serve as first-order approximations 
which ignore the uncertainties in the orbital parameters 
of eccentricity and argument of periastron and instead fo- 
cus on the time-domain parameters of period and time of 
periastron passage. However, the equations also serve to 
over-estimate the size of the transit window (the conser- 
vative approach) by assuming that the orbital inclination 
is edge-on compared with the line of sight. The conse- 
quence of this is that the maximum transit duration is 
allowed for. 

2.2. Orbital Parameter Dependencies 

As previously mentioned, the quality of the transit 
windows depends upon the uncertainties in the orbital 
parameters. The transit windows also grow with time, 



prompting follow-up of the transit window as soon as 
possible after discovery. In this section, we describe these 
effects for 245 known exoplanets for which the equations 
in Section 12.11 have been used to calculate their transit 
ephemerides. 

Figure [5] shows the size of the first transit window (the 
first transit to occur after t p ) and the uncertainty in the 
transit mid-point for the 245 exoplanets in the sample. 
Also shown are the predicted transit durations and geo- 
metric transit probabilities for these exoplanets. These 
calculations are based on the errors in the orbital pa- 
rameters derived from currently available radial velocity 
data. 

The necessity of the logarithmic scale in these plots 
demonstrates the large range in the size of the transit 
window. The transit windows of the short period plan- 
ets tend to be significantly smaller since, at the time of 
discovery, many orbits have been monitored to provide a 
robust estimate of the orbital period. This is particularly 
obvious in the plot of the transit mid-point uncertainics 
which shows that periods less than ~ 10 days have far 
superior constraints on the calculated ephemerides. In 
contrast, the longer period exoplanets often only have 
one orbit completely monitored and it is possible, though 
uncommon, for the resulting transit window to become 
comparable to the orbital period of the planet. The ideal 
targets to monitor in an observing campaign tend to oc- 
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Fig. 3. — The increase in the size of the transit window from the 
predicted time of the first window after t p to the predicted window 
after a JD of 2454979.5. The open circles shown in the bottom- 
right are those long-period planets for which an additional transit 
window beyond t p has not occurred as of this JD. 

cupy the lower-right corner of the plot of transit win- 
dow as a function of transit probability. These plan- 
ets have the highest likelihood of yielding successful de- 
tections, though this population is dominated by short- 
period planets. 

Figure [3] shows the net increase in the size of the transit 
windows for this sample of exoplanets by comparing the 
first transit windows after discovery with the first tran- 
sit window occurring after a JD of 2454979.5 (CE 2009 
May 28 00:00 UT). The open circles shown in the bot- 
tom right of the figure are those long-period planets for 
which an additional transit window beyond t p has not yet 
occurred and so the size of the transit window remains 
unchanged. Note that the distribution of points in this 
plot now resembles the distribution shown in the transit 
mid-point uncertainty plot of Figure [2J since the transit 
duration estimate is unaffected by the passage of time. 
Therefore, the transit window size increase for the short- 
period planets is much slower over time than for the long- 
period planets. This indicates that, even though many 
more orbits of the short-period planets have occurred, 
the transit mid-point uncertainty remains dominated by 
the uncertainty in the period. The size of the transit 
window for the long-period planets can be brought into 
a managable regime for photometric follow-up with rela- 
tively small usage of large telescope time. Without such 
an effort, it clear from these plots that it will be impossi- 
ble to ascertain whether or not many of the long-period 
planets transit their host stars. 

2.3. Improvements from Additional Data 

As described in Section 12.21 a considerable number of 
high transit probability targets are not feasible (depend- 
ing upon telescope access) to observe because the un- 
certainty in the predicted transit mid-point are too high 
to justify the observing time required. This can lead to 
transit windows of months and even years in duration. 
The acquisition of just a handful of new radial velocity 
measurements at carefully optimised times can reduce 
the size of a transit window by an order of magnitude. 
Here we describe, by way of two examples, how obtaining 
further radial velocity measurements for known exoplan- 
ets can improve the transit ephemerides. These examples 



were chosen based upon their very different periods, rela- 
tively high transit probabilities, availability of radial ve- 
locity data, and different transit windows and discovery 
dates. In each of the examples, we have simulated four 
additional measurements by using the best-fit orbital pa- 
rameters to determine the radial velocity at later epochs 
and adopting the mean of the discovery data precision 
for the simulated measurement uncertainties. The simu- 
lated measurements were then passed through a gaussian 
filter, which produced scatter consistent with the uncer- 
tainties, then appended to the discovery data. 

2.3.1. HD 190228 

The planet orbiting t he star HD 190228 was discovered 
bv lPerrier et~all ([2003) as part of a group of new planets 
announced by the ELODIE team. The planet is in a 
~ 1146 day orbit around a G sub-giant star with an 
eccentricity of ~ 0.5. The eccentric nature of the orbit 
resulted in no radial velocity data being acquired by the 
discovery team when the planet was close to periapsis, 
since the planet spends a very small portion of its orbit 
near that location. Calculations for the first predicted 
transit to occur after JD 2454979.5 (see Figure[3]) yield a 
transit mid-point uncertainty of 88.9 days and a transit 
window of 178.9 days. The geometric transit probability 
of this planet is ~ 1% which is relatively high for a planet 
of this orbital period. However, the large transit window 
makes this an unfeasible target to observe, particularly 
from the ground where a substantial fraction of the total 
transit window will remain uncovered (assuming only one 
ground-based telescope at a particular longitude is being 
used). 

In Figure |4] we show the discovery data of lPerrier et al.l 
(2003) along with four additional simulated measure- 
ments. The simulated data are each separated from each 
other by 50 days. Note that the simulated measurements 
have been acquired while the planet is speeding past pe- 
riapsis. The periastron passage of an orbit, particularly 
for a highly eccentric orbit, is where the planet is moving 
the fastest and so occupies a relatively small fraction of 
the total phase space. Thus, the greatest constraints dur- 
ing the shortest period of time can be made by sampling 
this part of the orbit. These effects have been discussed 
at length in the con text of the effect s of eccentric orbits 
on period analysis (Cumming 2004), cadence opt imiza- 
tion for radial velocity surveys dKane et al.l 120080 . and 
adaptive scheduling algorithms (|Fordll2008h . 

The orbital parameters of HD 190228 were re- 
computed from the combination of discove ry and sim- 
ulated dat a using the metho d described by iKane et al.l 
(|2007f) and IKane efaLl (|2009f ). The original and revised 
orbital parameters are shown in Table [TJ based upon the 
fits to the original and revised datasets respectively. The 
revised orbital parameters have only slight improvements 
in their uncertainties with the exception of the period 
and the time of periastron passage, which are the two 
most important parameters for calculating the transit 
window. This improvement from only four additional 
measurements decreases the uncertainty in the transit 
mid-point and the size of the transit window (for the 
first transit to occur after JD 2454979.5) by a factor of 
~ 6. Though the transit window is still quite large (31.0 
days), it is now far more accessable and of course can be 
improved further by increasing the phase coverage and 
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TABLE 1 
Fit parameters for HD 190228b. 



TABLE 2 
Fit parameters for HD 231701b. 



PetrfLmctcr 


Origin&l fit 


Revised fit 


P £iT ell I1C t C I' 


Origin&l fit 


Revised fit 


P (days) 


1146 ± 16 


1144.14 ± 2.09 


P (days) 


141.6 ±2.8 


141.89 ±0.15 


V (km s" 1 ) 


-50.182 ±0.004 


-50.181 ± 0.003 


V (m s" 1 ) 




-2.413 ± 1.824 


K (km s- 1 ) 


91 ±5 


90.75 ±4.36 


K (m s" 1 ) 


39.0 ±3.5 


39.06 ± 2.64 


u> (•) 


100.7±|| 


101.03 ±4.14 


w(°) 


46 ± 24 


54.40 ± 3.69 


e 


0-499±g:g|l 


0.501 ±0.041 


e 


0.10 ±0.08 


0.096 ± 0.069 


t p (JD - 2450000) 


1236 ± 25 


4672.076 ± 9.085 


t p (JD - 2450000) 


3180.0 ±4.2 


4885.141 ± 1.422 


td (days) 


1.155 


1.152 


t d (days) 


0.495 


0.491 


<5t mid (days) 


88.9 


14.9 


St mid (days) 


40.9 


1.6 


twin (days) 


178.9 


31.0 


twin (days) 


82.3 


3.7 



Note. — The original orbita l parameters for HD 190228b as 
measured by Pcrricr ct al. (2003) and the revised orbital param- 
eters with the four addditonal measurements, along with original 
and revised transit duration, transit mid-point uncertainty, and 
transit window. 



Note. — The original orbital parameters for HD 231701b as 
measured by Fischer ct al. (2007) and the revised orbital parame- 
ters with the four additonal measurements, along with original and 
revised transit duration, transit mid-point uncertainty, and transit 
window. 
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Fig. 4. — Best-fit solution (solid li ne) to the orig i nal ra dial veloc- 
ity data of HD 190228 obtained by IPerrier et all J2003T) and four 
subsequent simulated measurements. 

time baseline of the radial velocity data. The main bene- 
fit to constraining the transit window will come through 
improving the baseline (measurements during the same 
phase at subsequent orbits) rather than additional phase 
coverage, since phase coverage mainly aids toward con- 
straining the shape (eccentricity and periastron argu- 
ment) of the radial velocity variation. 

It should be noted that this analysis does not take 
into account the more typical situation where the addi- 
tional measurements are acquired with a different tele- 
scope and/or template spectrum than the discovery data. 
In this case, a floating offset (whereby the radial velocity 
offset between datasets is included as a free parameter) 
between datasets will need to be app lied during the fit- 
ting p rocess, such as that described by I Wright fc Howard! 
( 2009). This however has a neglible effect on the accu- 
racy of the fitted orbital parameters provided the addi- 
tional measurements have suitable phase coverage. 

2.3.2. HD 231701 

A more recent planet discovery is that of HD 231701b 
bv iFischer et al.1 ((2007) . This planet has an orbital pe- 
riod of ~ 141 days with an eccentricity of ~ 0.1. Even 
so, the slight eccentricity and an argument of periastron 
near 90° gives the planet an elevated geometric tran- 
sit probability of ~ 1.3%. The host star for this planet 
is a late-F dwarf. The data acquired at discovery was 
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Fig. 5. — Best- fit solution (solid li ne) to the orig i nal ra dial veloc- 
ity data of HD 231701 obtained bv IFischer et al.l (20071) and four 
subsequent simulated measurements. 

sufficient to constrain the orbital period to within a cou- 
ple of days. However, enough time has transpired since 
discovery such that the first predicted transit after JD 
2454979.5 has a mid-point uncertainty of 40.9 days and 
a total transit window of 82.3 days. 

Figure [5] s hows the discovery data published by 
IFischer et al.l (|2007f ) along with four additional simulated 
radial velocity measurements. The simulated data are 
each separated from each other by 10 days. As was the 
case for HD 190228, we found that the optimized con- 
straint on the period resulted from spacing the new mea- 
surements to cover a large range of radial velocity (am- 
plitude) space rather than phase space. This comes at 
the expense of refining the shape of the periodic varia- 
tion which, as described earlier, is determined by e and 
uj. The results from performing a fit to the combined 
dataset are shown in Table O The significant improve- 
ment to both the precision of the period and time of 
periastron passage parameters results in a subsequent 
improvement to the uncertainty in transit mid-point and 
transit window size that is impressive - a factor of almost 
25! This would result in the first transit window beyond 
JD 2454979.5 being a highly accessible window to obtain 
good coverage, particularly if longitude coverage could 
be achieved through appropriate collaborations. 

3. PHOTOMETRIC FOLLOW-UP STRATEGY 
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Amongst the southern hemisphere RV planets, there 
are few that have been adequately monitored photomet- 
rically to confirm or rule out planetary transits. By cal- 
culating transit ephemerides for these planets and design- 
ing an efficient observing program, it is possible to exam- 
ine a large subset of these planets with a relatively small 
amount of observing time on a l.Om-class telescope. By 
applying strict criteria on the predicted transit properties 
of the targets, we are able to produce a robust selection 
which yields the most promising targets on which one 
can place transit constraints. These criteria and general 
observing strategy considerations are described here. 

3.1. Target Selection 

For each observing run, the experimental design consti- 
tutes the selection of targets which meets the necessary 
criteria (described below) for successful observations to 
be undertaken. The primary challenge is to match tran- 
sit windows with observability for each target which is 
not a trivial task. Here we describe the minimum crite- 
ria that must be met for each target. The first step is to 
select all known radial velocity targets for which an es- 
timate of the stellar radius, either from measurement or 
models, is available for calculating the transit depth and 
the geometric transit probability. The probability that 
is most important for the target scheduling is the transit 
detection probability which is a combination of the geo- 
metric transit probability and the fraction of the transit 
window during which the target is observable. To first 
order, this is a straight multiplication but will depend 
upon the probability distribution of the predicted transit 
mid-point, as discussed in Section 13.21 Based upon the 
photometric precision of the experimental system, the 
predicted transit depth can be used to exclude targets 
whose depth is too low, particularly giant host stars. 

The steps thereafter depend upon how the observing 
time is allocated; a fixed time-slot (such as NOAO time), 
or queue-schedule d/ service time (such as that used by 
iKane et al.l ((2009)). For a fixed time-slot, the essential 
steps are: 

1. For 0.9-1.0m class telescopes, stars brighter than 
V = 6.0 often need to be excluded unless the tele- 
scope has the options of aperture diaphragms or 
neutral density filters. 

2. Include only those stars whose airmass is less than 
~ 2 for at least 3 hrs during the night. A transit 
window of 3 hours is the likely minimum transit 
window available and so the visibility of the target 
will still be useful if the transit window happens to 
largely coincide. 

3. For each target that passes the visibility and bright- 
ness criteria, the transit ephemeris is checked and 
transit times are noted for those predicted transits 
which fall on dates during the run. 

4. For each date during the run, the UT times of the 
transit are checked and transit windows which oc- 
cur during the day are rejected. Transit windows 
which do not coincide with the observable hours 
are also rejected. 

5. The schedule for each night is considered in terms 
of the transit detection probability (described 
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above) and the targets are ranked in descending 
order of the probability. If one seeks to concen- 
trate the investigation on long-period planets, then 
the planets should be ranked by the observability 
of the transit windows since the geometric tran- 
sit probability will dominate the transit detection 
probability for short-period planets. 

Queue-scheduled observations essentially allow an op- 
portunity to target long-period planets whose transit 
windows occur far less frequently than those usually 
monitored during pre-allocated observing runs. For 
queue-scheduled observations the steps are as follows: 

1. Rank the exoplanets from long-period to short- 
period and keep only those whose transit proba- 
bility exceeds the g eometric transit probabilit y for 
a circular orbit (see lKane fc von Brau nl lpOOl ). 

2. Examine the the transit window for each planet and 
reject those for which the transit windows are ex- 
cessively long (for example, greater than ~ 5 days). 

3. Investigate the visibility of each target from the ob- 
serving site and reject the targets for which there 
is a mis-match between the occurance of the tran- 
sit window and the time the target is observable. 
Ideal targets are those for which the transit win- 
dow occurs during a single night and the target is 
up all night. For targets which have transit win- 
dows spanning multiple nights, the decision will be 
based upon the value of the target in terms of the 
transit detection probability. 

These steps are almost the reverse of the steps recom- 
mended for fixed time-slot observing runs. Since queue- 
scheduled observations allow for targeted observations of 
rare transit windows for long-period planets, they are 
thus those which will yield the highest success for these 
high-risk/high-return targets. 

3.2. Coverage of Transit Window 

The difficulties in establishing an optimal observing 
schedule not withstanding, there are further considera- 
tions that one needs to take into account when planning 
observations. One of these is the decision about which 
parts of the transit window to monitor if the window 
spans more than one night (and often several nights for 
long-period planets), especially if one or more of those 
nights competes with other favourable targets. So far 
we have assumed that all parts of the transit window 
are equally significant in the likelihood of a transit be- 
ing observed or ruled out. This is generally not true but 
depends upon the probability distribution of the orbital 
parameter uncertainties. For example, we can suggest 
an empirical model whereby we assume Gaussian uncer- 
tainties for the fit parameters, which in turn assumes 
Gaussian noise in the radial velocity measurements. 

Figure [6] shows a Gaussian probability distribution 
for the predicted location of the transit mid-point for 
a planet which has a transit window size 5% of the or- 
bital period. If one is able to monitor the target for only 
half of the transit window, then choosing this range to 
be centered on the predicted transit mid-point (shown as 
the shaded region in Figure [6]) will account for 38% of 
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Fig. 6. — Probability distribution for the predi cted location of 
the transit mid-point (as described in Section I3.2H . assuming that 
the uncertainties in the orbital parameters are Gaussian. 

the area under the probability distribution, as opposed 
to 30% for the remainder of the time within the transit 
window. 

In reality, the uncertainties associated with the orbital 
parameters will have a more complex distribution due to 
systematic noise components. Their distributions may be 
close to that described by Gaussian or Poisson statistics, 
but can be determined empirically through Monte-Carlo 
simulations which randomize the sequence of the residu- 
als on the radial velo city measur emen ts and redetermin - 
ing the orbital fit (see lFordl (|2005D and lKane et~al] (|2007l ) 
for more details). In addition, the probability distribu- 
tion for a particular parameter is usually non-symmetric 
in nature, as seen in a \ 2 map of the parameter whilst 
keeping t he other pa r amete rs fixed (for example, see x 2 
maps by iKane et alj (|2009f) ). but can be approximated 
as symmetric at the la level. Whatever the distribution, 
the measurements per unit time will be more valuable the 
closer they are to the predicted mid-point. The caveat 
to this is when the predicted transit duration is larger 
than the time for which the target is observable during 
the night, since it is important to observe either ingress 
or egress for relative photometry. Thus, if the transit 
duration is greater than the observing window then the 
optimal approach is to observe at t m id ± td/2. 

The coverage of the transit window can be increased 
through the use of telescopes adequately separated in 
longitude. Such networks are already in existence (e.g., 
Las Cumbres Observatory Global Telescope (LCOGT) 
network) and collaborations for the follow-up of gamma- 
ray bursts and microlensing events are quite common. 
Queue scheduling of observations is particularly useful 
for the rare transit windows of long-period planets. This 
kind of observing is available, for instance, to member 
consortiums who utilise the service time of the Observing 
with Small and Moderate Aperture Research Telescope 
System (SMARTS). 

For the especially bright targets which will not only 
saturate typical l.Om-class telescopes but whose field-of- 
view will also be devoid of comparison stars, a solution 
is to use the Micr ovariability and Osc i llations of Stars 
(MO ST) satellite (jWalker et all I2003t iMatthews et all 
I2004D . The MOST satellite has demonstrated photomet- 
ric precision of a few parts per million which is suffi- 
cient for detecting transit signatures due to planets orbit- 



ing bright giant stars. Since MOST is space-based, this 
would also allow complete coverage of the transit window 
without the need for coordinated ground-based observa- 
tions using different telescopes (an additional source of 
red noise) . Provided the transit window can be provided 
with sufficient accuracy, this would be an excellent use 
of the MOST satellite's capabilities. 

4. CONCLUSIONS 

Many of the known radial velocity planets have yet to 
be surveyed for transit signatures. The detection of a 
transit for the intermediate to long-period planets would 
add enormously to our knowledge of planetary struc- 
ture and, in particular, how the structure varies with 
semi-major axis and periastron distance. The advan- 
tages of targeting long-period radial velocity planets are 
the brightness of the host stars and the prior knowledge 
of the planetary orbital parameters. However, a major 
challenge of monitoring the host stars at predicted tran- 
sit times is that many transit windows have deteriorated 
over time, such that the telescope time required renders 
attempts to do so impractical. 

We have shown through calculations for 245 of the 
known exoplanets how the size of the transit window 
varies with period and geometric transit probability. The 
large uncertainties associated with the transit mid-point 
for the long-period planets is dominated by the uncer- 
tainties in the period and time of periastron passage esti- 
mated from the discovery data. We demonstrated, using 
the examples of HD 190228 and HD 231701, that a hand- 
ful of carefully timed additional measurements can vastly 
improve the size of the transit window and thus bring the 
monitoring of the window into the reach of ground-based 
programs. 

The difficulties involved in the observing schedule 
largely result from matching transit windows with the 
observability of the targets, particularly for long-period 
planets whose transit windows are widely spaced. We 
have described a planning strategy which will make opti- 
mal use of both pre-allocated and service telescope time, 
also noting the advantages of both longitude coverage 
and space-based observations. It is important to con- 
sider for the scheduling that the central part of the tran- 
sit window can be significantly more valuable than the 
wings of the window, depending upon the nature of the 
orbital parameter uncertainties. 

The described techniques and science goals are 
currently being undertaken and investigated by the 
Transit Ephemeris Refinement and Monitoring Survey 
(TERMS). Note that the observations from this survey 
will lead to improved exoplanet orbital parameters and 
ephemerides even without an eventual transit detection 
for a particular planet. The results from this survey will 
provide a complimentary dataset to the fainter magn i- 
tude range of the Kepler mission (|Borucki et alj [20091. 
which is expected to discover many transiting planets 
including those of intermediate to long-period planets. 
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